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Abstract 

An exact solution of the two center time-dependent Dirac equation for pair 
production induced by ultrarelativistic heavy ion collisions is presented. Cross 
sections to specific final states approach those of perturbation theory. Multi- 
plicity rates are reduced from perturbation theory. 
PACS: 25.75.-q, 34.90. +q 

I. INTRODUCTION 



In this paper, we shall compute the production of electron-positron pairs in the central 
region for highly relativistic charged ions. We base our calculations on an exact solution 
of the time dependent Dirac equation in the ultrarelativistic limit. We show that except 
for possible cutoff effects the exact cross section for any specific final electron positron 
state equals the perturbation theory result. On the other hand, we argue that the rate for 
processes correlated with low impact parameters (such as pair multiplicity) are reduced from 
the perturbation theory rate. 
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Our notation will vary at times between that of high energy physics (light cone variables 
and Dirac 7 matrices) and that of atomic physic (Dirac a. and f3 matricies and the usual 
four momentum and space time variables) as convenience and the connection to previous 
work dictates. The relation between different notations should be clear. 

To properly define what we mean by central region and highly relativistic, it is useful to 
define light cone coordinates, 

P ± = ^( Po ± Pz ), x ± = ^=(t±z) (1) 

where the z direction is the beam direction. In this coordinate system, the invariant dot 
product of momenta is p ■ q = p + q~ + p~q + — px • qr- We will take the nucleus propagating 
along the positive z axis to have a large light cone momentum P + and that along the negative 
z axis to have large light cone momenta Q~ . By highly relativistic we mean, for example, 



the case of RHIC, colliding ions each with 7 (= 1/vl — v 2 ) of 100. Note that the Compton 
wavelength of the electron is large compared to the radius of either nucleus and thus the 
nucleus can be considered a point charge aZ. 

To define what we mean by centrally produced, we define the lightcone momentum 
fractions of an electron or positron to be 

x = p+/p+ (2) 

and 

y = P~/Q- (3) 

Central production will mean electrons and positrons which have x,y « 1 This can be 
satisfied for some range of longitudinal momenta of the electron or positron so long as the 
condition on the Compton wavelength is satisfied. 

When the electron or positron in the central region sees the moving nuclei, it sees two 
oppositely Lorentz boosted Coulomb fields. These are the Lienard-Wiechert potentials 
= aZ(l-va z ) + aZ(l + va z ) 

v/[(b/2- P )hY + {z-vty ^/[(b/2 + P )hY + {z + vtf 
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b is the impact parameter, perpendicular to the z-axis along which the ions travel, p, z, 
and t are the coordinates of the potential relative to a fixed target (or ion) a z is the Dirac 
matrix, and Z, v and 7 are the charge, velocity and 7 factor of the oppositely moving ions. 
We have specialized to the case of equal Z ions; unequal Z would only require a trivial 
change in what follows. If one makes the gauge transformation on the wave function 

V> = e- ix(r 'V (5) 

where 

®Z , r , . rz ~ ~ aZ 
— In 7(2 - vt) + Jb 2 + -f 2 (z - vt) 2 ] 

V v V 



X(r, t) = — ln[7(z - vt) + ^ 'b 2 + 7 2 {z - vt) 2 } - — ki[j(z + vt) + ^b 2 + 7 2 (z + vt) 2 ] (6) 



the interaction V(p, z,t) is gauge transformed to |] 

aZ(l — va z ) aZ(l — (l/v)a z ) 



V(p,z,t)- 

^[(b-pv^ + ^-vt) 2 v /w + 5 - 

aZ(l + va z ) aZ(l + (l/v)a z ) 



y/[(b + p)h} 2 + (z + vt) 2 v /6 2 / 7 2 + (z + ^) 2 ' 

This gauge transformation reduces the range in (z ± vt) to more closely map the B and 
E fields (which have the denomenator to the | power rather than the \ power of the 
untransformed Lorentz gauge). In the ultrarelativistic limit (ignoring correction terms in 
[(b + p)hf) Eq.(7) takes the form § 

V(p, z, t) = 5(z -t)(l- a z )A-(p) + 5(z + t)(l + a z )A + (p) (8) 

where 

A ± (p) = -Z„h6^, (9) 

The potential as written here will be referred to the singular gauge solution. In this gauge 
the field vanishes everywhere except along the lightcone x ± = 0. We can gauge transform 
from this field to the less singular light cone gauge by again utilizing ip = e~ lx(j,t ^il)' where 

X (r, t) = 9{t - z)A-(p) + 6(t + z)A + (p). (10) 



This leads to added gauge terms in the transformed potential 



dx(r,t) 
dt 



- a ■ V 



6{z-t)(l 



a z )A-(p) - S(z + *)(! + a z )A + (p) 



6{t-z) et 



VA-(p)-6(t + z) a- VA + (p) 



(11) 



and we thus obtain the light cone gauge 



V(p, z,t) 



6{t -z)a- VA-(p) - 6{t + z) a ■ VA+(p) 



(12) 



with a. the Dirac matrix. This construction fits in well with the similar non-abelian treat- 



The reason why we expect that we can exactly compute pair production in the central 
region should be obvious. In either of the gauges above, the propagation of the electron 
or positron is essentially trivial. Except at x = 0, the electron or positron propagates 
as a free particle. Therefore, to construct the propagator for the electron, one needs to 
solve a boundary value problem with free propagation everywhere except at the surfaces of 
discontinuity at the lightcone x ± = 0. 

In fact to solve the problem of pairs production is arguably a little simpler than con- 
structing the propagator. What we will do is assume we begin with a negative energy 
solution of the Dirac equation in the initial state. This will correspond to a positive energy 
positron. We will then let it propagate forward in time. At late times we will compute the 
amplitude that this state is a positive energy electron. 

The organization of this paper is as follows: In the second section, we compute the 
amplitude for electron positron pair production in the singular (5 function) gauge. The 
technique used will closely follow that previously used in the exact calculation of bound 
electron positron pair production and will use conventional Green's functions methods. 
In the third section we work in light cone gauge. We show that the result for the pair 
production amplitude agrees between the two computations. In the fourth section, we 
discuss the cross sections and their relation to perturbation theory. In Appendix A, we 
write down our conventions for light cone coordinates and projection operators in the Dirac 



ment for QCD §. 
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equation. In Appendix B, we evaluate the transverse integral of the Coulomb field which is 
necessary for the solution in both gauges. In Appendix C, we show how to solve the Dirac 
equation across a boundary corresponding to a charged nucleus in light cone gauge. We 
show how this result maps into the similar result computed in singular gauge. 

II. SOLUTION IN SINGULAR GAUGE FOR PAIR PRODUCTION 

A strategy to find the exact semi-classical solution for pair production in a two-center 
model is to first find the exact wave function at the point of interaction (in terms of the 
appropriate Green's function) and then to construct the exact amplitude (incorporating the 
initial exact wave function, the interaction, and the final plane wave function). 

In particular, one would like to solve the two center time-dependent Dirac equation, 




(13) 



in terms of a Dirac Green's function based on the plane wave solutions, 




(14) 



One obtains 



op + P - E (r, t) = V(p, z, f)tf (r, t)+ ap + (3 — E 0(r, t) 



(15) 



If one acts on both sides with the Dirac plane wave Green's function 



i -l 



(16) 



G = ap + [3 - E 



one obtains the usual form 



q? = ( f) + G Vy. 



(17) 



The Green's function G obeys the equation 



ap + (3 - E G q (t, t; r', t') = S(r - r')S(t - t') 




(18) 
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and the solution is 

-y p gjfc(r— r')— ibj(t— t') 

Go(r, t- r', = (2^1 y [«k + /3 + w ] fc2 - i _ ^ . (19) 

Note the positive sign of uj in the brackets. We will get the boundary conditions later. 

We now make use of the usual light cone coordinates, x + , x~ , q + , and q~, as in Eq.(l) in 
place of t, z, go, and q z . Begin with a plane wave representing one of the electrons in the 
filled negative energy Dirac sea. Our convention will be that go is positive, corresponding 
to the positron energy we eventually are interested in, and likewise the three momentum is 
opposite in sign to that of the negative electron 

= v ( q> Sl )e- iqx = v(q, Si ) e -^^ q +x-+i q -x+ _ ^ 

The function v(q, s$) is the usual spinor of a positron. The Green's function can likewise be 
written in light cone coordinates 

G (x; x') = ^4 J dk + dk-d 2 k ± [oik + (3 + u] _ 2 k+k~ + k\ + l ' (21) 

and in light cone coordinates the two center time-dependent potential Eq.(8) takes the form 
V(x) = -^=5(a:-)(l - a z )A-(x ± ) + ^(x + )(l + a z )A+(x ± ). (22) 

A space-time diagram is presented in Fig. 1. Obviously V(x) only acts on the boundaries 
between the Regions I, II, III, and IV. 

In Region I we have the initial plane wave 0. We would like to begin by constructing 
the solutions in Regions II and III. We will work out Region III explicitly and then Region 
II will follow by symmetry. We have from Eq. (17) 

*m{x) = 4>(x) + J d 4 zG {x, z)V(z)*{z). (23) 

But we have on the boundary between I and III @] 

V(z)V(z) = -4(500(1 - a z )A-(z ± )e i6 ^ A ~^(t)(z), (24) 
V2 



which can be expressed equivalently 



V(zMz) = ^S(z-)(1 - a z )(e iA ~^ - !)4>{z). 



(25) 



We now obtain 



^iii(x) = 4>{x) - J dz + dz d 2 z±j^j J dk + dk d 2 k ± 



x 



(2tt)< 
1 1 

a±k± + f3 + — ^=(1 - a z )k + —j=(l + a z )k~ 



ik±(xj_—zj_)—ik + (x —z )—ik (:r+— z+) 

-2k+k~ + k\ + 1 



x -l=S(z-)(l - a z ) (e lA ~ (2i) - Si ) e - iq ^ +iq+z ~ +iq ~ z+ . 

The term in (1 + a z ) obviously vanishes. Now integrate over z~, z + , and k~ to obtain 



(26) 



dk+d 2 k , 



a±k± + (5 Wl — a z )q~ 

v2 



D ik±x± — ik+x +iq x^ 

2k+q~ + k\ + 1 



x-^=(l-a z ) | rf 2 Zie -(«i+'i)^( e 'A-(^) _ i) v ( 5> s .). 



(27) 



Now consider the integration over k + . We did not previously specify the boundary 
condition on Go and we must do so now. We want the Go term in Eq.(27) to be non- 
vanishing for x~ > and to vanish for x~ < 0. If we move the singularity just below the 
real axis by adding ie to the denominator then this physical boundary condition will be 
fulfilled. We have 



-c 



dk+e -ik+x 

k+ + ^— + ie 



-2irie i[(k ^ +1)/(2q ~ )]x ~ 



(28) 



and our expression for the wave function in Region III now becomes 
f d 2 k 1 

^iii(x) = - J [-7=(a±k ± +P)(1- a z ) - (1 - a z )q~ 



x ■ 



(2n)^V2 

e ik ± x ± +i[(k 2 ± + l)/(2q-)]x-+iq-x- t 



/ 



c p z±e -i(l±+k±)zj_ e iA (z ± ) 



v(q, s^. 



(29) 



As a check on the validity of our expression we note that if we project with (1 — a z ) the 
expected boundary condition holds for x~ = +e: 



{l-a z )^ HI = (l-a z )e lA '^U 



(30) 



Now the exact amplitude takes the form 

M = i J dt < <f) \V \ty > . (31) 

Begin by constructing V times the exact wave function on the boundary between Regions 
III and IV. In analogy to Eq. (25) we have 

V(x)V(x) = ^5{x + ){l + a,)(e lA+(lJ - l)* m (a:). (32) 
v2 

Note that we need the projection (1 + a z ) for this boundary due to the opposite direction 
of motion of the ion producing the V as compared to the V on the boundary between I and 
III. Thus we have at x + = 0, x~ > 

V(xMx) = ^(x+)(l + a z y^ J H^K^ + m - «,) 

ik ± x ± +i[(k 2 L +l)/(2q-)]x-+iq-x+ f 

x- — J d 2 z ± e- i{ ^ +k ^e ih ~ iz ^v{q, Si ). (33) 

Likewise at x~ = 0, x + > 

V(x)*(x) = j=8(x-)(l - a z )e iA ~^ J ^^( a± k ± +(3)(1 + «,) 

ik ± x ± +i[{k 2 ± + l)/(2q+)]x++iq+x- f 

x- — d 2 z ± e- i ^ +k ^e iK+ ^ ) v{q 1 s i ). (34) 

Zi(J J 

In constructing the transition amplitude one makes use of the fact that only two in- 
teraction (two photon) terms have a net contribution. (Single interaction terms integrated 
over the four boundaries give a null contribution). The amplitude then has two pieces 
corresponding to the boundary of Region IV with Regions II and III 

M(p, q) = M(p, q) IVJII + M(p, q) IVJI . (35) 

The final state is a positive energy electron 

= u (p, s f )e ipx = u{p, Sf yp±x±-ip + x--ip-x + . (36) 

Then we have 



roo r 

M(p,q)ivjn — i J dx~ J dx + dxi_u(p, Sf)/3< 



, 3 ip±x±+ip + x +ip X -1 



ik ± x ± +i[(k 2 ± +l)/(2g _ )]x _ +ig _ x + f 

x — / d 2 z ± e~ i ^ +k ^e iA ^v{q, Si ). (37) 

Now integrate over x + and x~ to obtain 

d 2 k± u(p, Sf)f3(l + «z)(a^A;^ + f3)v(q, Sj) 



m(p, q)iv,iii = - J 



(2tt) 2 2p+q~ + k\ + 1 

x y rf 2 a; ± e- i(px - fex):c± e iA+(:c±) ^ dVe- i(<?x+fex) *V A ~ (z±) . (38) 

The transverse spatial integrals can be done in closed form as is shown in Appendix B. We 
obtain 

J y± ~ W k 2 ~ 2iv r(iv) [ j 



where 77 = +Za. 

The amplitude now becomes 



M(p,, Wr =-4-i--M/ 



r 2 (l — i?7) u(p, s/)/3(l + a«)(o!±fcL + /3)u(g, 



T 2 (ir]) J 2p+q~ + k 2 ± + l 

Xe -i6(px/2- 9 x/2-fex)(( p± _ A; ± ) 2 (A; ± + g^) 2 )^ 1 (40) 

where trivial b dependent and constant phases have been removed. 

The other piece of the amplitude coming from the boundary between Region II and 
Region IV has the corresponding form 

T 2 (l-iri) r 2 tx(p, s/)/?(l - a z )(a±k ± + (3)v(q, Si ) 
M(p, q) IVJI = -4 v2{ir]) J d k ± 2p _ g+ + kl + 1 

Xe +»(px/2-,x/2-fcx)((p_ L _ A; ± ) 2 (A; ± + gx) 2 ) 1 "" 1 
1 r2 ( 1 ~^) /"^ - a*)(a±(p± - 5± - k±) +(3)v(q,Si) 



T 2 (i7]) J ' 2p q+ + (p ± - q ± — k ± ) 2 + 1 

xe- lb(p ^ 2 -^ /2 -^\( Pl _ - k ± ) 2 (k ± + q±) 2 ) ir >-\ (41) 



and the total amplitude takes the form 
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M(p,q)=4 V 2 J d 2 k ± ((p ± - k±) 2 (k± + q ± f ! '"' 

'u(p, Sf)P(l + a z )(a±k± + ff)v(q, si) 



2p+q~ + k]_ + 1 

u(p, s f) /3(1 - a z )(a±{p± -g±- k±) + /3)v(q, sj) \ , 

2p-q+ + (p ± - g ± - k ± f + 1 j ' " 



where trivial phases depending on r\ and on initial and final momenta have been removed. 

Rewriting with the mass of the electron explicit, and making use of more modern notation 
we have 



M(p, q) = Arf J d 2 k ± e ibk ± ((p ± - k ± ) 2 (k ± + qA _ff ' 
'u(p,s f )(l- a z )(- ^ ± + m)v(q,Si) 



2p + q~ + k 2 _ + m 2 
+ M(p,s/)(l + a z )(- d±+ ft± + m)v(q,Si) " 
2p~q + + (p_i_ — gj_ — A;_|_) 2 + m 2 



At this point one notices that the infinite range of the transverse potential provides 
us with a result lacking an infrared cutoff. In Appendix B, we introduce a convergence 
parameter e to regularize this infrared singular behavior. Strictly speaking we were not 
allowed to let the convergence parameter e take on the value zero. In the previous use of the 
singular interaction for calculating bound electron positron pair production, the transverse 
integrals were cut off by the finite range of the bound state wave function [|J . With plane 
waves there is no such natural cutoff. This suggests that we attach a physical interpretation 
to the convergence factor e in Appendix B. In fact if we set e equal to uj/j, where u is the 
energy of the produced electron or positron, then we obtain the expected spatial cutoff in 
a heavy ion electromagnetic interaction, and a result that can naturally be compared with 
the corresponding perturbation theory result of Bottcher and Strayer 0. We discuss this 
further in Section IV. But first we shall come to the same solution in the light cone gauge. 
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III. THE SOLUTION TO THE TWO CHARGE PROBLEM IN LIGHT CONE 

GAUGE 



In this section, we present an alternative derivation to the results presented above in 
singular gauge. This serves as a check on the results of the previous section. Also the 
technique which we present here appears to generalize directly to the non-abelian problem. 

We begin by considering a negative energy plane wave. Here we will take q ± < 0. This 
negative energy incoming plane wave will eventually be interpreted as an incoming positron 
with momentum q — > —q. So we begin with an incoming wave 

Mx) = e tqx v x (-q). (44) 

Here V\ is a positron spinor. The label I on the solution refers to the fact that this is a 
solution in region / as in Fig. 1, and as in Section II. 

Recall that the background field in light cone gauge is of the form 

A\x) = -6(x-)V l A-{x T ) - 6{x + )V l A + {x T ) (45) 

as was discussed in the introduction. 

Let us first construct the solution in region 77. This involves constructing the solution 
in region across the boundary at x + = from region I. Everywhere in region 77, x~ < 0. 
Using the results of Appendix C, a general solution is of the form 



lh TT (x) = _L e -i^ + M f d P +d2 PT c {-ip+x-+ip T x T -i(pl+m*)x+ /2(p++ie)} L , a T ' PT + m 

V2 J (2tt) 3 \ y/2(p+ + ie). 



a + v x (-q) (-— l - / d 2 z T e l(qt ~ PT)zT e tA+(zT) + H(p T ,p + ,q)\ (46) 

[i p + — q + — ze J ' J 

In this equation, H has all its singularities in the negative half p + plane. Therefore at x + = 0, 
this contribution gives nothing. In fact throughout region 77, there is no contribution so 
that H can be dropped. The choice of singularity for p + in the exponential is so that there 
is convergence at large positive x + . Notice that for x~ < corresponding to region II, the 
term which does not involve H may be evaluated by closing in the upper half p + plane. 
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Here the solution satisfies the correct boundary condition at x + = and solves the Dirac 
equation in region II. It is unique. 

In exactly the same way, we have in region that the solution is 

rh TTT (x) = —e-^-M f d P~ d2 PT c {-ip-x++i PT x T -i( P l+ m i)x- /2( P -+ie)} L , «t ; Pt + m \ 
V2 J (2vr) 3 \ y/2(p- + ie) J 

oTv x {-q) {-— - / d 2 z T e l{qt - pT)zT e lA ~ M ) (47) 



i p~ — q~ — ie 

We now write down an ansatz for the solution in region IV. After writing down the 
solution, we will show that it solves the Dirac equation and the boundary conditions at 
x ± = when x T > 0. Consider 

^J X ) - i e -iA+(x T )-iA-( XT ) f d 2 PT d 2 Z T d 2 p' T d 2 z' T 



dp + 1 



2n p + — + m 2 )/2q- + ie 

^~ip+x~+ip' T x T —i(p' 2 +m 2 )x + /2(p++ie) 
e iz' T (p T -p' T )+iz T (q T -p T ) e *A+ (z T )+iA~ (z' T ) 



' a T -p' T + pm\ ,a T -pT + l3m _ 

1 H =- — -= a V\(—q) 

V2(p+ + ie) J y/2q~ 



, I dP \ e -ip- x+ +ip' T x T -i(p'2+m 2 )x+ /2(p- +ie) 

J 2n p~ - {p\ + m 2 ) /2q+ + ie 



e iz' T (p T -p' T )+iz T (q T -p T ) e iA {z T ) +i A+ (z' T ) 

First it is straightforward to apply the Dirac equation to ipiv, and show that it solves 
the Dirac equation. What about the behavior on the boundaries x ± = 0, x T > 0. First 
consider x~ = and x + > 0. We first evaluate ipn on this boundary. We require that the 
ip~ be continuous. Using ip n , we see that 



aT ' P ^ + (3m a + v x (-q) ( d 2 z T e^-^e* A+ ^ (49) 

Now, how does vfjjy behave on this boundary? The first integral over p + in the integral 
representation for ipj v vanishes at x~ = as the all the singularities of the integration over 
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p + are in the same side of the real p + axis in the complex p + plane. In the second term 
which involves p~ , we must close in the lower half p~ planes when x~ = and x + > 0. 
Closing in this plane and doing the integrals over z' T and p' T shows that our formula solves 
the correct boundary condition. 

By symmetry, we see that along the other boundary at x + = for x~ > the boundary 
condition is solved for the plus component of the wavefunction. Therefore iprv solves the 
Dirac equation and the boundary conditions in region IV. 

Now we must extract the matrix element for pair production from this wavefunction. To 
do this, we observe that after a little algebra, and making the substitution q — > —q that 



WJ " 2 6 J (2nY (2n)* d ZTd Zt 

m — p' • 7 



-e~ tp x | 
p' 2 + m 2 — ie 



2n' + 

e iz' T (p T -p' T )-iz T (q T +p T ) 



P l+ + (Pt + m 2 )/2q~ + ie 
e iA +W ^- Ma -rn + PT-7 Vx{q) 



V2q~ 



p 



2p' 



iz l T {p T -p' T )-iz T (q T +p T ) 



p'- + (Pt + m 2 )/2q+ + ie 

e *-W+* + M a + m + PT-7 Vx ( q) } (50) 
\ / 2q + p' 

(51) 



Amputating the external line which corresponds to the propagator for the electron, we 
find the matrix element for pair production to be 



M " X = ljik { K+iP ' T " PTh ' iP ' T + qT) 2q-^ + vl + m 2 + i e U ^ a ' {m + P ' T ' 7K(g) 



(2tt) 2 l '" xrr xri ' ^'2q-p'+ +p' 2 + m 2 + ie 

+n + (p' T + q T )n-{p' T - p T ) - * - ^ - - u^p) V2a + (m + p' T ■ j)v x (q) } (52) 

Here the quantity 

«±(p r ) = J d 2 x T e- ipTXT e +lA±(xT \ (53) 
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IV. DISCUSSION 



After we had finished the derivation of the amplitudes Eq.(43) and its equivalent Eq.(52) 
we became aware of similar results recently obtained by Segev and Wells pi. These authors 
obtained an amplitude with the same structure as ours, but they did not obtain our closed 
form, Eq. (40) (B6) or Eq.(B5) for the non-perturbative transverse momentum integral. 
They did argue that the transverse integral went appropriately to the ultrarelativistic per- 
turbative limit in agreement with the result of Bottcher and Strayer ||. In this limit their 
perturbative amplitude simply takes the form of Eq.(43) without the ir] in the exponent. 

We now have the suggestion of an ansatz for a physical cutoff of the exact result in line 
with the discussion at the end of Section II. The perturbation theory result of Bottcher 
and Strayer not taken to the ultrarelativistic limit actually has denominators of the form 
kj_ + u 2 /(j 2 — 1) rather than k\. Our ansatz then is to modify Eq.(43) to 

M(p, q) = 4 V 2 [ d 2 k ± e«*± ((\p± - k ± } 2 + ^) (> ± + q ± } 2 + ^) 

, f u{p,s f ){l - a z ){- ji± + m)v(q, Sj) 
V 2p + q~ + k 2 L + m 2 

u(p, s f ) (1 + a z ) (- 4 L + ft± + m)v(q, Sj) \ ^ 
2p~q + + (p± — q± — k±) 2 + m 2 J 

We have effectively retained e (set tou/j) in the factor taken to the 1 — irj power of Eq.(B5) 
and ignored it elsewhere. Retaining rj elsewhere in Eq.(B5) would cut off a little more 
sharply at small k± but is not necessary. It would not necessarily be more exact to retain 
all factors of rj because the cutoff comes in response to the spatial region p = j/u where 
both the singular and light cone potentials begin to lose their validity. 

At this point one is left with a choice whether to perform the integral over the impact 
parameter either before or after the integral over the intermediate transverse momentum. 
If one is only interested in a the cross section for pair production in a given small momen- 
tum and/or energy bin then one must obtain the same answer independent of the order 
of integration. In that case performing the integral over the impact parameter first seems 
convenient. One has an impact parameter integral of the form 
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J d 2 h I \M(p,q)\ 2 ~ Jd 2 hJ d 2 k ± e* bk -f(k ± ) J d 2 k' ± e- ibk '±r(k' ± ). (55) 

The integral over d 2 h gives (27r) 2 <5^ 2 ^ (A:^ — k' ± ) and the ±ir) exponents vanish, giving a result 
identical to what we would obtain in perturbation theory. 

If we are interested in high multiplicity events or a total cross section for events in 
which at least one pair is produced then we must perform the integral over the intermediate 
transverse momentum first. For the highest multiplicity events we need to evaluate the 
square of the amplitude at small but non-intersecting values of the impact parameter: we 
need the mean number of pairs at a given impact parameter. From the mean number of 
pairs produced at each impact parameter, one can obtain the probability (less than one) of 
at least one pair being produced. Thus a total pair production cross section can be defined 
and computed. 

Note that for a given impact parameter the square of the amplitude is not identical in the 
exact solution to what it is in the the perturbation theory solution. For the exact solution 
the % rj in the exponent gives a rotating phase in the the integral over k± that is absent in 
the perturbative case. For example if we look at a dominant contribution near one of the 
cutoffs the amplitude goes as 

M ~ I* — k±d 9 k± ^ . = J-^^WVt 2 ) _ W/t 2 ) (56) 
Jo (kl + iu 2 /^ 2 ) 1 -^ 2irj y 1 v 1 

which decreases with increasing 77. Here the exact contribution is less than perturbation 
theory would be. 

To sum up, we have a situation where one can say both that the exact cross section for 
pair production to any final state is identical to the perturbation theory cross section, and 
that measurements of high multiplicity pair events will show deviations from perturbation 
theory probably by being smaller. 
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APPENDIX A: CONVENTIONS FOR LIGHT CONE VARIABLES AND THE 

USE OF PROJECTION OPERATORS 

We will let the light cone coordinates for any vector be 

V ± = ±{V°±V). (Al) 

The metric convention we will use is g 00 = 1 and g' lj = —5 tj for spatial components so that 
in light cone gauge the dot product is 

A ■ B = A+B- + A-B + - A T ■ B T . (A2) 

We shall use ordinary gamma matrices for fermions so that 

{y,y} = 2«r. (A3) 

We can define projection operators for fermions as 

P ± = ^l°l ± = 1 -(l±a z ). (A4) 

These projection operators satisfy P + + P~ = 1, P ±2 = P ± , and P ± P T = 0. The plus and 
minus components of a fermion field are 
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^± = (A5) 

where ip + + ip~ — 

In light cone coordinates, the Dirac equation for a free particle becomes 

^ = (pr ■ oct + Pm)^ + (AG) 

where %p + solve 

(2p + p- -p T - m 2 )tp + = 0. (A7) 

In an external field which is independent of x , the transverse momentum is converted 
into a transverse covariant momentum. Also, the relationships can be reversed under the 
transformation ip — > 

APPENDIX B: EVALUATION OF THE TRANSVERSE INTEGRAL 

We would like to evaluate the integral 

d 2 p e- lk - p e v wv 2 . (Bl) 



A coordinate shift gives us 

(I? \ * n r2n roc 

— J e Ti[k - b/2] J d(j) J pdp e -(ifc±cos^+ e |cos0|)p e -i n lnp ? ( B2 ) 

where we have added a convergence factor e\ cos0|. The integral over p can then be carried 
out ( H 3.381.5) and we obtain 

b2 \ ir ' „ T i[k-b/2] f 2n . . r ( 2 ~ (— 2i— 2rj) arctan (fcx cos </>/ej cos </>[) 



e Ti[k-b/2\ / V U c (-li-Zr)) arctan (fc x cos tp/e\ cos 0|) (B3) 

4 7 io [(A;i + e 2 )cos 2 0] 1 - i " ' v ; 

Making use of quadrant symmetry this may be rewritten 



I = 4 e 



k . b / 2 ]/6 2 y' 7 (e 2 — A; 2 ) cosh(2?7 arctan [fcx/e]) + «2e/cj_ sinh(2r/ arctan [/cj_/e]) 
Tj (e 2 + k 2 ) 



r(2 - 2in) W 2 d<j> 

(k 2 ± + e 2 ) 1_ "> Jo (cos0) 2 - 2 ^' 1 ' 
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The integral over can then be carried out by Wallis' Formula ( 6.1.49) giving 

=F»[k-b/2] (b 2 \ ir > (e 2 — k\) cosh(2r] arctan [kj_/e\) + i2ek± sinh(2r] arctan [k±/e\) 

U) (e 2 + kj) 

r(2 - 2»77) 0Fr(-.5 + ^) 

(fci + e 2 ) 1 "^ 2r(itj) ' 1 J 

Letting e approach zero from the positive direction and exploiting T function relations finally 
leads to 

UeJ ^- 2i " r(<f/) ' 1 j 



APPENDIX C: SOLVING THE DIRAC EQUATION FOR THE 
LIENARD-WIECHERT POTENTIAL 

In this appendix, we discuss solving the Dirac equation for Lienard-Wiechert potentials 
appropriate for the central region. We will show how to construct a solution for a light cone 
potential of the form 

A l T = -^(x")V i A(x T ) (CI) 

We will then show how to convert this solution into that for the singular gauge and how this 
problem in singular gauge translates into a boundary condition on the Dirac wavefunction 
at x~ = 0. In the text, we will use this to construct solutions across the boundaries at 
x ± = 0. 

We assume that for x~ < 0, we are given a plane wave solution of the form 

i/>(x) = e iqx u x (x) (C2) 

where A is a polarization label on the electron spinor. At x~ > 0, the solution must be a 
linear combination of plane wave solutions to the Dirac equation modulo a gauge rotation 

^'( x ) = e * A ^) J 0LF xx ,(p T , q T )u' x (p)e^ XT -^ x+ -^ + ^ (C3) 



where 

P + = ^ (C4) 

The integral above is two dimensional since q~ labels the solution on both sides of the 
boundary as the potential is x + independant, and because the value of plus component of 
light cone momenta is determined by the mass shell condition. 

What are the boundary conditions at x~ = 0? Recall that the form of the Dirac equation 

is 

(P T -aT + pm- V2p + P- - ^/2p-p + )^ = (C5) 

where P ± are the projection operators described in the last appendix. Here Pt is the 
covariant momentum operator. We see there for the ip~ can be chosen to be continuous 
across the boundary. On the other hand, the plus component of the wavefunction must be 
discontinuous since by the Dirac equation 

i) + = -^—{P T -a T + l3m)i)- (C6) 

and the covariant momentum operator contains the vector potential which is discontinuous. 

This tells us what the proper boundary condition is in singular gauge. Recall that to 
transform between the two gauges 

1^ singular{%) ^ ^ ^ T ' ^ ^lightcone(%) ■ (^-''^) 

We therefore have that the boundary condition in singular gauge is the the component 

^vngularix) U"=0+= ^ H'singular^) U"=0- (C8) 

with the plus components determined by the Dirac equation. All components of the wave- 
function in singular gauge are therefore discontinuous. 

It is now straightforward to determine the solution for x" > 0. Consider 
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= e * A ( g T) / ^ ^ T d 2 Zj>c~ iA ( ZT " > c ipTXT ~ iq x+ ~ i P +x e iz T(QT-PT) 



V2 J (2tt^ 2 

•Pt - 
V2q 



a T -pT + l3m \ _ 
H t=— fa « A (g) (C9) 



where 



a ± = 4=(l±«z) = v / 2P ± . (CIO) 
v2 

Again p + = (pf, + m 2 )/2q~ It is easy to see that this function solves the Dirac equation. At 
x~ = 0, the term involving p + disappears so that the integrals over px and zt can be done 
with the result 

= ^ (a" + a + ar -^!' 8m ) e * U »( 9 ). (Oil) 

Using the definition of the projection operators P 1 * 1 , and the relationship between this 
is just e iqx u\(q). To derive this, we must use the definition of the vector potential in terms 
of A. This solution therefore solves the boundary conditions. 
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FIGURES 

FIG. 1. Light cone boundaries of the four regions of the z, t plane for an ultrarelativistic 
collision. 
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